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Abstract 

Based on results of density functional theory (DFT) calculations with the local 

ri ' spin density approximation (LSDA) and the generalized gradient approximation 

—L \ (GGA), we propose a new magnetic material, CsCl-type FeSe. The calculations 

C , reveal the existence of ferromagnetic (FM) and antiferromagnetic (AFM) states 

O . over a wide range of lattice constants. At 3.12 A in the GGA, the equilibrium 

state is found to be AFM with a local Fe magnetic moment of ±2.69 ur- A 
K* ' metastable FM state with Fe and Se local magnetic moments of 2.00 ub and 

(N : 

— 0.032 /XB, respectively, lies 171.7meV above the AFM state. Its equilibrium 
iy-\ . lattice constant is ~ 2 % smaller than that of the AFM state, implying that 

lO ' when the system undergoes a phase transition from the AFM state to the FM 

o 

f~^ I one, the transition is accompanied by volume contraction. Such an AFM-FM 

transition is attributed to spin-density z-reflection symmetry; the symmetry 
driven AFM-FM transition is not altered by spin-orbit coupling. The relative 
5_^ ' stability of different magnetic phases is discussed in terms of the local density 

of states. We find that CsCl-type FeSe is mechanically stable, but the magnetic 
states are expected to be brittle. 
Keywords: CsCl-type FeSe, first-order magnetic transition. 
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1. Introduction 

Iron-selenide shows a number of polymorphs including a low-temperature 
tetragonal (a) phase, the so-called anti-PbO-type FeSe p|. The a-Fe^Sei-^ 
compound is currently an interesting phase and has been studied experimentally 
for its spintronics related magnetic properties [2, y, |j, l5| . The perfectly stoichio- 
metric a-FeSe has been shown to be non-ferromagnetic, but non-stoichiometric 
a-FeSe demonstrates ferromagnetism that was attributed to defects or Fe clus- 
ters in a-FeSe thin films |5|. Zinc blende- type FeSe was shown to be an anti- 
ferromagnetic metal [6| . Recently, superconductivity with a Tc '^ 8 K was also 
discovered in a- FeSe and was attributed to Se vacancies I7f|. Subedi et al. j8| 
performed first-principles calculations and their results indicated that electron- 
phonon coupling can not explain the superconductivity of a-FeSe which shows 
spin-density wave (SDW). Zhang et al. [9] showed that excess Fe in FeTe pro- 
vides electron and the excess Fe is strongly magnetic. 

Motivated not only by the observed magnetism and superconductivity in 
a-FeSe, but also by our recent extensive calculations of FeSe 111], where we con- 
sidered different crystal structures of FeSe and (i) a phase transition from a-FeSe 
to CsCl-type FeSe was found and (ii) CsCl-type FeSe was the only structure that 
survived in a compressed unit cell and retained its magnetism. The tetragonal 
a-FeSe structure with coordination number four can be thought of as the CsCl 
structure with an elongated c axis and coordination number of eight. Since 
ordered CsCl-type systems have no nearest neighbor atoms of the same kind, in 
order to understand the properties of such compounds containing magnetically 
active constituents, it is necessary to consider magnetic unit cells large enough 
to allow for antiferromagnetism. 

Here, we investigate the magnetism of CsCl-type FeSe. We find that CsCl- 



type FeSe demonstrates both fcrromagnctism and antiferromagnetism over a 
wide range of lattice constants in terms of density-functional theory (DFT) [lO| 
by using the total-energy all-electron full-potential linearized augmented plane- 



package 



m 



wave (FLAPW) [ll| calculations implemented in the QMD-FLAPW software 



2. Computational Method 



Both the local spin density approximation (LSDA) [13| and the generalized 
gradient approximation (GGA) [IJl were considered for the exchange-correlation 
potentials. Integrations inside the Brillouin zone (BZ) were performed using 
the improved tetrahedron method [15j over a 15 x 15 x 15 mesh within the 
three dimensional (3D) BZ, corresponding to 120 k points inside the irreducible 
wedge of the 3D-BZ. An energy cutoff at 4.1231 (27r/a), where a is the lattice 
constant, was employed for the linearized augmented plane- wave (LAPW) basis 
set, which corresponds to ^ 210LAPWs per k-point and spin. A 16.3707 (27r/a) 
star function cutoff was used for depicting the charge density and potential in the 
interstitial regions. Lattice harmonics with I < 8 were employed to expand the 
charge density, potential, and wave- functions inside each muffin-tin (MT) sphere 
of radius 2.2 a. u. for Fe and 1.9 a. u. for Se. Note that those computational 
parameters satisfy the convergence test [la |. 

All core electrons were treated fully relativistically and valence states were 
calculated scalar relativistically, i.e., without spin-orbit coupling [17| . The ex- 
plicit orthogonalization (XO) scheme was employed to ensure the orthogonality 
between the core and valence states [l8|. For spin-orbit coupling on valence 
states, we employed the second variation method [19| with the spin diagonal 
parts of the density subjected to a self-consistency loop. During the second 
variation procedure, integrations inside the 3D-BZ were done in the fuU-BZ, i.e., 
with 1688 k points. Self-consistency was assumed when the difference between 
input and output charge densities became less than 1.0 x 10^^ electrons/a. u.'^ 
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Figure 1: The antifcrromagnctic unit cells for (a) type-I, (b) type-II, (c) type-Ill, and (d) 
type-IV presented on the basic CsCl unit cell for FeSe. Open and filled circles represent Fe 
and Se atoms, respectively. Arrows indicate different spin directions of the Fe atoms. 



3. Results and Discussions 



3.1. Magnetic Phase Transition 

Because the LSDA and GGA results agree qualitatively with each other, 
we will keep our discussion based on the GGA results; the LSDA results will 
be given where they are required. The total energies were calculated for the 
nonmagnetic (NM) and ferromagnetic (FM) states using the conventional unit 
cell of two atoms. Figure [1] shows the antifcrromagnctic (AFM) coupling of Fe 
in different planes and four types of AFM structures, AFM-I, AFM-II, AFM- 
III, and AFM-IV, were considered and tetragonally doubled unit cells were used 
for the AFM calculations. The main difference among these AFM structures is 
not only the AFM coupling between the Fe atoms in different planes, but also 
antiferromagnetically coupled coordination numbcr(j/) in the unit cell e.g., the 
AFM-I, AFM-II, AFM-III, and AFM-IV states have two, four, six, and four 
antiferromagnetically coupled atoms, respectively. Furthermore, AFM-II and 
AFM-IV have the same coordination number, but the coupling between the 
Fe atoms is different when viewd in the xy plane. Note that all these AFM 

n 

structures are similar to the FeAs-based superconductor [20|. 

The calculated total energy energy curve shown in Fig. [2] demonstrates that 
the AFM-I state is the ground state of CsCl-type FeSe. The calculated total 
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Figure 2: The calculated total energy in units of inRy/Fe atom as a function of lattice constant 
(A) of FeSe in the NM, FM, AFM-I, AFM-II, AFM-III, and AFM-IV states relative to the 
equilibrium total energy of the AFM-I state. Filled circles, squares, and triangles present the 
NM, FM, and AFM-I states, respectively. Open circles, squares, and triangles represent the 
AFM-II, AFM-III, and AFM-IV states of CsCl-type FeSe, respectively 



energies per Fe atom are shown relative to the equilibrium energy of the AFM-I 
state. The total energy difference between the FM and AFM states is sensitive to 
the unit ceU volume change. Since the AFM-II, AFM-III, and AFM-IV states are 
higher in energy than the FM and AFM-I states, we denote for brevity hereafter 
the AFM-I state as the AFM state. The calculated equilibrium lattice constants 
for the NM, FM, and AFM states are summarized in Table[TJ Calculations were 
also carried out for the optimization of the c/a ratio of the AFM structure, and 
the optimized c/a was found to be ^ 1.985. 



Table 1: Calculated properties of CsCl-type FeSe; a is the equilibrium lattice constant in A 
units, mpo is the local magnetic moment of Fe in fi^ units, and AE gives the relative stability 
of the FM, AFM-II, and NM phases with respect to the AFM-I phase equilibrium in units 
of meV/Fe atom, i.e., AE < means that the AFM-I phase is more stable than the FM, 
AFM-II, and NM phases. Results are given for LSDA and GGA. 





LSDA 


GGA 


Phase 


a (A) 


mpo im) ^E 


(meV/Fe) 


a (A) 


rriFc (p-b) 


A^ (meV/Fe) 


NM 


2.94 


- 


-229.97 


3.01 




-480.59 


FM 


2.96 


1.70 


-40.78 


3.05 


2.00 


-171.70 


AFM-I 


3.00 


2.12 


0.00 


3.12 


2.69 


0.00 


AFM-II 


2.96 


1.37 


-192.55 


3.05 


2.09 


-341.38 



The energy difference per Fe atom between the AFM and FM states, 
^^ ^ i;;(AFM) - i;(FM) 



(1) 



where n is number of Fe atoms, reflects the inter-atomic exchange coupling 
strength between Fe spins at the AFM equilibrium lattice, is calculated to be 
about -205 meV (-56meV) with GGA (LSDA). The AFM state is more stable 
than the FM state for lattice constants larger than the critical value, ac = 2.90 A, 
where the AFM-FM transition occurs. It is to be mentioned that the LSDA 
gave flc = 2.92 A. The magnetic state of CsCl-type FeSe is very sensitive in the 
region 2.90-2.94 A and this is the region where FeSe makes a transition from the 
AFM state to the FM one under compression. The energy barrier {^E) can be 
overcome if we decrease the lattice constant by a ~ 7 % (~ 2.6 %)( °°''~°^ x 100) 
with GGA (LSDA) to resuh in the AFM CsCl-type FeSe in the FM state, as 
seen in Fig. [2j To be more confident, we used the same tetragonal unit cell and 
found the FM stability against AFM in the region where the FM is more stable 
than the AFM. 

It will be more interesting if we know how much external pressure is necessary 
to achieve this metastable CsCl-type FeSe. As we found, CsCl-type FeSe has 
different magnetic structures and one can calculate the transition pressures of 
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Figure 3: Calculated local Fe magnetic moments (in /^b units) within each MT sphere. The 
inset shows the the local magnetic moments of Se for FM and AFM states. Squares and circles 
represent the local magnetic moment of Fe for the FM and AFM states, respectively. 



these magnetic structures. However, we only calculated the transition pressures 
for the NM, FM, and AFM states. The transition pressure was calculated 
by using the equation of state [2l| which enables us to calculate the enthalpy 
H = E + pV , where E is the internal energy, p is the external pressure, and V 
is the volume of the system. The transition pressure from a-FeSe to CsCl-type 
FeSe is obtained from the usual condition of equal enthalpies, i.e., the pressure 
at which the enthalpies H of a-FeSe and CsCl-type FeSe are the same. We 
followed this approach and the transition pressure is determined by equating 
the enthalpies of the two phases of FeSe. The GGA shows that pressure values 
of about 24.32, 20.45, and 19.90 GPa are necessary for the transition of a-FeSe 
to NM, FM, and AFM CsCl-type FeSe, respectively. 

The calculated local magnetic moments within each MT sphere of the Fe 



and Se atoms are shown in Fig. [3] and the local Fe magnetic moments of all 
the magnetic phases are summarized in Table [1] For both the FM and AFM 
states, the local magnetic moments of Fe increase as one increases the lattice 
constant. It is noticeable that there is a crossing point around 2.79 A. Below 
this lattice constant, where the FM state is more stable than the AFM one, 
the Fe local magnetic moments in the FM states are larger than those of the 
AFM state. Above 2.79 A, the Fe atom has a larger magnetic moment in the 
AFM state than in the FM state. The induced Se magnetic moment in the FM 
state was found to be -0.032 /^b (-0.012 ^b) with GGA (LSDA). In Fig.[3l we 
can see that in the FM region the induced magnetic moments at the Se site are 
positive, whereas in the region where the AFM state is the most stable state, 
Se has negatively induced magnetic moments. A slight jump of the Fe moments 
in the FM state can be seen near the critical region of 2.90-2.94 A, where the 
transition takes place. These local magnetic moments show that the magnetic 
state (FM or AFM) is stabilized by the high spin state of the Fe atom. 

Fe Se Fe Se Fe 




Figure 4: Spin-density contours in the (110) plane of CsCl-type FeSe for the (a) FM and 
(b) the AFM states. The solid and broken lines represent spin-up and spin-down densities, 
respectively. The lowest contour starts from 2 X 10~* electrons/a. u.'' and the subsequent lines 
differ by a factor of V2. The Fe and Se atoms are located at the corner and center of each 
panel, respectively. 



It is interesting to find tliat the calculated local magnetic moment of the Se 
atom in the AFM state is zero. To understand this feature, the spin-density 
contour plots of CsCl-type FeSe are shown in Fig. Ufor both the FM and AFM 
states. The FM state shows, for the spin-density z-reflection operation, even 
symmetry, while the AFM state shows odd symmetry. This implies that the 
FM-AFM transition should be a first-order phase transition, because a first- 
order phase transition is understood as the choice of the representation of a 
symmetry operation 



22|. Kim et al. [23| revealed the spin-density- inversion 
symmetry driven first-order FM-AFM phase transition of GaCMua. Considering 
the fact that the inversion operation can be understood as a rotation followed 
by a refiection, it is possible to generalize this idea to the spin-density-refiection 
symmetry driven first-order FM-AFM phase transition. 

In order to check the effects of SOC, we also carried out calculations in the 
FM and AFM states with SOC. However, there is no significant alteration of 
the symmetry driven transition by SOC. 

3.2. Electronic Structures 

The relative stability among the magnetic orderings can be understood from 
an analysis of the calculated density of states (DOS). Figure [SJa) shows the 
DOS for the NM state of CsCl-type FeSe; the Fe d states are decomposed into 
t2g and Cg states. Note that the high peak at the Fermi level {E-p) originates 
from Fe-Cg states hybridized with the Se-p states. The large DOS at E-p, which 
comes mainly from Fe-e^ states, indicates the ferromagnetic instability in terms 
of the mean-field Stoner theory [2J|. One may observe that the t2g states lie 
between the Cg states. Such behavior is due to the structure of CsCl-type FeSe 
which is typical for body-centered cubic metals where a pseudogap separates 
the bonding and anitbonding states [25|. Such a splitting can also be seen in 
the FM DOS. Once the FM ordering happens, as in Fig. [5fb), one can see the 
exchange splitting in the Fe-d states, along with the fact that the Fe-egf states 
are occupied completely, while the Fe-t2gt states are partially occupied. Strong 
hybridizations between the Fe-d states and the Se-p states are also observed in 
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Figure 5: Total and atom-projected density of states of the (a) NM, (b) FM, and (c) AFM 
states of CsCl-type FeSe at their equihbriuni lattice constants. Solid, dotted, and dashed lines 
show Fe t2g, Eg and Se p states, respectively, whereas solid lines in the bottom panels show 
the total DOS. The Se p states are multiplied by a factor of 10 and the total density of states 
of AFM FeSe is given per formula unit for comparison purposes. The Fermi (Ep) levels are 
set to zero. 



both spins. 

An important difference between the FM and AFM structure is that in FM 
the effective field split the bands and hence the peaks in the total DOS, while 
in the AFM case as seen in Fig. [SJc) , the spin- up and spin-down peaks do not 
exhibit any such splitting and differ only in their relative intensity. An additional 
feature of the AFM phase is the absence of Se local magnetic moments consistent 
with the spin density maps. The populations of the partial states of Se atoms 
are very small in both magnetic structures, therefore the Fe d states which are 
dominant at Ep and behave differently in the FM and AFM phases are primarily 
responsible for their relative stability. Furthermore, in the AFM structure, Ep 
does not fall in the valley, the majority Fe-t2g and Fe-e^ states are completely 
occupied, and the DOS at Ep is very small and leads to a stable AFM structure, 
while the FM splitting leads to an increased DOS at Ep mainly contributed by 
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Figure 6: Total energy per Fe atom (in units of eV) versus (a) tetragonal and (b) orthorhombic 
strain of FeSe in their NM, FM, and AFM phases. The total energy is given with respect to 
unstrained AFM FeSe for tetragonal and orthorhombic distortion. Filled solids, squares, and 
triangles represent the NM, FM, and AFM-I states, respectively. 



Fe Fe-i2g- In addition, the Fe-d spin- up states are mostly occupied and provide 
the relatively higher local Fe magnetic moment to be 2.69 /iB, which is larger 
than that of the FM case of 2.00 /xb. state. It is to be noted that AE, mFe, and 
density of states at Ep of CsCl-type FeSe are larger than those in a- FeSe [8] . 

3.3. Mechanical Stability 

Finally, an important question to address is whether this new material is 
mechanically stable or not. To answer this question we calculated the tetrag- 
onal and trigonal shear constants, C" and C44, respectively of the NM, FM, 
and AFM states, by computing the change in energy of the FeSe phase under 
small volume-conserving strains. The results of tetragonal and orthorhombic 
strains are shown in Fig. [B] and the calculated values are summarized in Table [5J 
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Table 2: The calculated bulk modulus {B), shear moduli (C and C44), and shear moduli on 
the slip plane (fi) in GPa units of CsCl-type FeSe in the NM, FM, and AFM (=AFM-I) states 
for both LSDA and GGA. 





LSDA 


GGA 


Phase 


B 


C 


C44 


A^ 


^i/B 


B 


C 


C44 


M 


/i/S 


NM 


174.6 


74.7 


73.8 


74.2 


0.42 


142.6 


66.6 


61.4 


63.5 


0.45 


FM 


146.6 


99.4 


74.0 


84.2 


0.57 


93.0 


80.1 


55.2 


65.2 


0.70 


AFM 


117.0 


90.3 


82.5 


85.6 


0.73 


82.7 


60.3 


61.7 


61.1 


0.74 



Considering the fact that CsCl-type FeSe is metastable [l|, it is mechanically 
stable in the NM, FM, and AFM phases against tetrahedral and orthorhombic 
deformations. The mechanical stability is confirmed by the positive value of 
the cubic elastic constant C" — (Cn — Ci2)/2 shown in Table [2j However, one 
should care on the brittle nature of the magnetic states as represented by the 
shear modulus [26| on the slip plane /i — (3C44 + Cn — C12) /5, to the bulk 
modulus B in Table [21 where all the magnetic states are predicted to be brittle 
^/B > 0.5. Note that the GGA predicts higher brittleness of the FM state than 
the LSDA; this fact reflects the difference in bonding nature of the approxima- 
tions. It is interesting to find that the NM state is also mechanically stable and 
rather ductile, i.e., 11/ B < 0.4, in contrast to bcc Fe [27|. It is to be noted that 
the detailed mechanical behavior of CsCl-type FeSe at high pressure (reduced 
lattice parameters) is beyond the scope of the current study. At reduced lattice 
parameters, the nature of the brittleness of CsCl-type FeSe may change. 

The lattice constant over which the FM state becomes stable is close to that 
of FM bcc Fe, 2.87 A, which has been grown on GaAs 



281 



29| andZnSe |30|,l31| 



substrates due to the good lattice match between them. Therefore, it can also 
raise the possibility of growing CsCl-type FeSe on GaAs or ZnSe, because the 
calculated lattice parameter of FeSe is very close to that of GaAs (oafm — 
aGaAs/2). The formation energy of CsCl-type FeSe at ambient pressure is ~ 
leV/f.u. HI, which indicates that CsCl-type FeSe is in a metastable state and 
there might be a possibility to achieve this structure either at high pressure 
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or grown as a thin film. Recently, FeS,which is isoelectronic to FeSe, indeed 



showed CsCl structure at high pressure 
have been grown successfully [33|, |3J, [35 , 



3211 . Several other nictastable materials 



36, 



37| using molecular beam epitaxy 



(MBE), which is a non-equilibrium process. 

4. Summary 

In summary, we predict, based on density functional theory by using the 
total-energy all-electron FLAPW method, a new phase of FeSe, the FM and 
AFM state of CsCl-type FeSe . The calculated total energy curves indicate that 
the AFM state is most stable with the largest lattice constant. In addition, the 
energy barrier for the AFM to FM transition can be overcome by volume con- 
traction. We found a spin-density-reflection symmetry driven first-order phase 
transition in CsCl-type FeSe, as seen previously in the spin-density-inversion 



symmetry driven phase transition in GaCMna [23[. The calculated LDOS and 
spin-density contour plots revealed the origin of the magnetic phase transitions 
in CsCl-type FeSe. We also predict that CsCl-type FeSe is mechanically stable 
for all the magnetic states considered in this work. 

Acknowledgements 

The authors appreciate H. K. D. H. Bhadeshia for his careful reading of 
the manuscript. I. G. Kim thanks Jong-Hoon Chung for helpful discussions on 
experimental realization. This work was supported in part by the Steel Innova- 
tion Program of POSCO, by the Basic Science Research Program (Grant No. 
2009-0088216) through the National Research Foundation funded by Ministry 
of Education, Science and Technology of the Republic of Korea, and by the U.S. 
Department of Energy (Grant No. DE-FG02-88ER 45372). 



13 



References 

[1] G. Rahman and I. G. Kim, Bull. Korean Phys. Soc. (Autumn Meeting) 26 
(2008) 66. 

[2] Q. J. Feng, J. Y. Zhang, B. S. Li, B. H. Li, Y. M. Liu, X. W. Fan, and H. 
W. Liang, Appl. Phys. Lett. 88 (2006) 012505. 

[3] K. W. Liu, J. Y. Zhang, D. Z. Shen, C. X. Shan, B. H. Li, Y. M. Lu, and 
X. W. Fan, Appl. Phys. Lett. 90 (2007) 262503. 

[4] X. J. Wu, D. Z. Shen, Z. Z. Zhang, J. Y. Zhang, K. W. Liu, B. H. Li, Y. 
M. Liu, B. Yau, D. X. Zhao, B. S. Li, C. X. Shan, X. W. Fan, H. J. Liu, 
and C. L. Yang, Appl. Phys. Lett. 90 (2007) 112105. 

[5] X. J. Wu, Z. Z. Zhang, J. Y. Zhang, B. H. Li, Z. G. Ju, Y. M. Liu, B. S. 
Li, and D. Z. Shen, J. Appl. Phys. 103 (2008) 113501. 

[6] G. Rahman, S. Cho, and S. C. Hong, J. Korean Phys. Soc. 53 (2008) 380. 

[7] F. -C. Hsu, J. -Y. Luc, K. -W. Yeh, T. -K. Chen, T. -W. Huang, P. M. Wu, 
Y. -C. Lee, Y. -L. Huang, Y. -Y. Chu, D. -C. Yan, and M. -K. Wu, Proc. 
Nat'l. Acad. Sci. U.S.A. 105 (2008)14262. 

[8] A. Subedi, L. Zhang, D. J. Singh, and M. H. Du, Phys. Rev. B 78 (2008) 
134514. 

[9] L. Zhang, D. J. Singh, and M. H. Du, Phys. Rev. B 79 (2009) 012506. 

[10] P. Hogenberg and W. Kohn, Phys. Rev. 136 (1964) B864; W. Kohn and L. 
J. Sham, Phys. Rev. 140 (1965) A1133. 

[11] E. Wimmer, H. Krakauer, M. Weinert, and A. J. Freeman, Phys. Rev. B 
24 (1981)864, and references therein; M. Weinert, E. Wimmer, and A. J. 
Freeman, tbid. 26 (1982) 4571. 



[12] For details see http://www.flapw.com 



14 



[13] U. von Barth and L. Hedin, J. Phys. C: Solid St. Phys. 5 (1972) 1629. 

[14] J. P. Pcrdew, K. Burke, and M. Ernzcrhof, Phys. Rev. Lett. 77 (1996) 3865; 
ibid. 78 (1997)1396(E). 

[15] J. -H. Lee, T. Shishidou, and A. J. Freeman, Phys. Rev. B 66 (2002) 233102. 

[16] S. -W. Seo, Y. Y. Song, G. Rahman, L G. Kim, M. Weinert, and A. J. 
Freeman, J. Magnetics 14 (2009) 137. 

[17] D. D. Koehing and B. N. Harmon, J. Phys. C 10 (1977) 3107. 

[18] M. Weinert, G. Schneider, R. Podloucky, and J. Rcdingcr, J. Phys.: Con- 
dens. Matter 21 (2009) 084210. 

[19] R. Wu, D. Wang, and A. J. Freeman, J. Magn. Magn. Mater. 132 (1994) 
103. 

[20] E. Aktiirk and S. Ciraci, Phys. Rev. B 79 (2009) 184523. 

[21] F. Birch, Phys. Rev. Lett. 71 (1947) 809; F. D. Murnaghan, Proc. Nat'l. 
Acad. Sci. U.S.A. 30 (1944) 244. 

[22] L. D. Landau and E. M. Lifshitz, Statistical Physics 3rd. Ed. (Elsevier, 
Amsterdam, 2007). 

[23] I. G. Kim, Y. J. Jin, J. I. Lee, and A. J. Freeman, Phys. Rev. B 67 (2003) 
060407(R). 

[24] D. A. Papaconstantopoulos, W. E. Pickett, B. M. Klein, and L. L. Boyer, 
Phys. Rev. B 31 (1985) 752. 

[25] J. Kiibler, Theory of Itinerant Electron Magnetism (Oxford, 2000). 

[26] Y. Wu and W. Hu, Eur. Phys. J. B 60 (2007) 75. 

[27] H. C. Herper, E. Hoffmann, and P. Entel, Phys. Rev. B 60 (1999) 3839. 

[28] J. S. Claydon, Y. B. Xu, M. Tselepi, J. A. C. Bland, and G. V. Laan, Phys. 
Rev. Lett. 93 (2004) 037206. 



15 



[29] J. J. Krcbs, B. T. Jonkcr, and G. A. Prinz, J. Appl. Phys. 61 (1987)2596. 

[30] M. Marangolo, F. Gustavsson, M. Eddrief, Ph. Sainctavit, V. H. Etgens, 
V. Cros, F. PctrofF, J. M. George, P. Bencok, and N. B. Brookes, Phys. 
Rev. Lett. 88 (2002) 217202. 

[31] J. J. Krebs, B. T. Jonkcr, and G. A. Prinz, J. AppL Phys. 61 (1987) 3744. 

[32] N. Sta, H. Ohfuji, K. Hirose, H. K. Ashi, Y. Ohishi, and N. Hirao, American 
Mineralogist, 93 (2008) 492. 

[33] H. Akinaga, T. Manago, and M. Shirai, Jpn. J. AppL Phys. 39 (2000) 
L1118. 

[34] J. F. Bi, J. H. Zhao, J. J. Deng, Y. H. Zheng, S. S. Li, X. G. Wu, and Q. 
J. Jia, Appl. Phys. Lett. 88 (2006) 142509. 

[35] J. H. Zhao, F. Matsukura, K. Takamura, E. Abe, D. Chiba, and H. Ohno, 
Appl. Phys. Lett. 79 (2001) 2776. 

[36] C. S. Tian, D. Qian, D. Wu, R. H. He, Y. Z. Wu, W. X. Tang, L. F. Yin, Y. 
S. Shi, G. S. Dong, X. F. Jin, X. M. Jiang, F. Q. Liu, H. J. Qian, K. Sun, 
L. M. Wang, G. Rossi, Z. Q. Qiu, and J. Shi, Phys. Rev. Lett. 94 (2005) 
137210. 

[37] M. G. Sreenivasan, K. L. Teo, M. B. A. Jahl, T. Liew, T. C. Chong, and 
A. Y. Du, IEEE Trans. Magn. 42 (2006) 10. 



16 



